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The velocity dependence of the stopping power of swift protons and 
deuterons in low energy collisions is investigated. At low projectile ener- 
gies the stopping is mainly due to nuclear stopping and charge exchange of 
the electron. The second mechanism dominates after E p > 200 eV. A dy- 
namical treatment of the charge exchange mechanism based on two-center 
electronic wavefunctions yields very transparent results for the exchange 
probability. We predict that the stopping cross sections vary approxi- 
mately as Vp 35 for projectile protons on hydrogen targets in the 1 keV 
energy region. 
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Nuclear fusion reactions proceed in stars at extremely low energies, e.g., of the 
order of 10 keV in our sun |1], @. At such low energies it is extremely difficult 
to measure the cross sections for charged particles at laboratory conditions due to 
the large Coulomb barrier. One often uses a theoretical model to extrapolate the 
experimental data to the low-energy region. Such extrapolations are sometimes far 
from reliable, due to unknown features of the low-energy region. E.g., there might 
exist unknown resonances along the extrapolation, or even some simple effect which 
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one was not aware of before. One of these effects is the laboratory atomic screening 
of fusion reactions || ||. It is well known that the laboratory measurements of low 
energy fusion reactions are strongly influenced by the presence of the atomic electrons. 
This effect has to be corrected for in order to relate the fusion cross sections measured 
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in the laboratory with those at the stellar environment. Another screening effect, 
arising from free electrons in the stellar plasma, will not be treated here. For about 
one decade, until 1996, one observed a large discrepancy between the experimental 
data and the best models available to treat the screening effect. The simplest (and 
perhaps the best of these models), the so-called adiabatic model, predicts that as the 
projectile nucleus penetrates the electronic cloud of the target the electrons become 
more bound and the projectile energy increases by energy conservation. Since the 
fusion cross sections increase strongly with the projectile's energy, this tiny amount 
of energy gain (of order of 10-100 eV) leads to a large effect on the measured cross 
sections. However, in order to explain the experimental data, it is necessary an extra- 
amount of energy - about twice the value obtained by the adiabatic model. This is 
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puzzling, since more refined dynamical models, e.g., time-dependent Hartree-Fock ||, 
include electronic excitation and thus yield a screening energy which is smaller than 
that obtained with the adiabatic model. 

This problem was apparently solved in 1996 by Langanke and collaborators || 
and by Bang and collaborators |7[|, who observed that the experimental data for 
3 He(d, p) A He - the reaction for which the screening effect was best studied - was 
probably obtained with a erroneous extrapolation of the stopping power for deuterons 
in helium targets to the low energy regime. The fusion reaction occurs at a point 
inside the target after the projectile has slowed down by interactions with the atomic 
targets. In the experimental analysis one needs to correct for this energy loss in 
order to assign the right projectile energy value for that reaction. These corrections 
were usually based on the Andersen-Ziegler table of stopping power of low energy 
particles [|J. Due to the lack of experimental information on the stopping power at 
the extreme low projectile energies needed for astrophysical purposes, the Anderson- 
Ziegler tabulation was extrapolated to the required energy; another example of a 
dangerous extrapolation procedure. In fact, Golser and Semrad || observed a strong 
departure of their experimental data from the extrapolations based on the Andersen- 
Ziegler tables for the stopping of low energy protons on helium targets. Grande 
and Schwietz [10] performed a dynamical calculation of the energy dependence of 
the stopping power for this system and confirmed that the extrapolation procedure 
cannot be extended to the very low energies. Whereas at higher energies the stopping 
is mainly due to the ionization of the target electrons, at the astrophysical energies 
it is mainly due to charge-exchange between the target and the projectile. Refs. || 
and use these arguments to explain the long standing discrepancy between theory 
and experiment for the low energy dependence of the reaction 3 He(d, p) A He. Other 
reactions of astrophysical interest (e.g., those listed in by Rolfs and collaborators 
0. M) should also be corrected for this effect. 

In this work we address the problem of the stopping of very low energy ions in 
matter. To simplify matters, we study the system p+H, which is the simplest one can 
think of. It displays important features of the stopping power and has the advantage 
of allowing a very simple solution. 

Our approach is based on the solution of the time- dependent Schrodinger equation 
for the electron in a dynamical two-center field. The static two-center p + H system 
has been solved by Edward Teller in 1930 ]TT[. He showed that as the distance between 
the protons decreases the hydrogen orbitals split into two or more orbitals, depending 
on its degeneracy in the two-center system. Analogous problems are well known in 



quantum systems ||12| . For example, take two identical potential wells at a certain 
distance. For large distances the states in one well are degenerated with the states 
in the other potential well. As they approach this degeneracy is lifted due to the 
influence of barrier tunneling. Thus, the lowest energy state of hydrogen, Is, splits 
into the lscr and the 2pcr states as the protons approach each other. The lscr state 
is space symmetrical, while the 2pcr state is antisymmetric. As the proton separation 
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distance decreases their respective energies decrease. At R ~ 1 A the energy of the 
2pcr state starts to increase again, while the energy of the lscr state continues to 
decrease. For proton distances much smaller than 1 A the lscr and the 2p<r energies 
correspond to those of the first and second states of the He atom, respectively [[11] . 



Let us now consider the dynamical case. The full time-dependent wavefunction 
for the system can be expanded in terms of two-center states, <fr n (t), governed by the 
Schrodinger's equation 

[•Ho + V p (t)] <t> n {t) = E n (t) <t> n (t) , with H = p 2 e /2m e + V T , (1) 

where V p (t) = — e 2 / |r + R/2| is the electron-projectile proton interaction potential 
and Vt = — e 2 / |r — R/2| is the electron-target proton interaction for a proton-proton 
separation distance R(t). Note that in our formulation the two-center wave functions 
also depend on time, as well as their energies E n (t). The full electronic wavefunction 
is obtained by a sum over all orthonormal two-center states 

\® (*)) = a « (*) (*)) ' with / d3r ^ (*) ^ (*) = Snm ■ ( 2 ) 

n J 

Inserting this expansion in eq. @j we obtain 

ifi i a m (*) = E m (t) «m (*) -ih^a n (t) ( m \-^\ n ) ■ ( 3 ) 

n 

Using (|I|) one can easily show that, for m ^ n, 

i i d , . (m\dV v /dt\n) . . , . 

(ml — \n) = — - — — , (m ^ n). 4 

v 1 dt l ' E n (t)-E m (ty K r } K J 

Moreover, using the second relation of eq. (0), one can show that (m\ -| \m) = 0, 
if \m) is real. This indeed will be our case. Our basis, \n(t)), is formed by two- 
center states at a given time t, i.e., a given proton separation distance, R. These 
wavefunctions are real. Thus, the final coupled-channels equation for the two-center 
problem is given by 

d , . , . v-^ (ml dVp/dt In) 

i/i— a m (t) = E m (t) a m (t) - ^ 2^ a « (*) E ( t )_E (t) ' ^ 

At very low proton energies (E p < IkeV) it is fair to assume that only the low-lying 
states are involved in the electronic dynamics. Only at proton energies of order of 25 
keV the proton velocity is comparable to the electron velocity, v e ~ ac. Thus, the 
evolution of the system is almost adiabatic at E p < 10 keV. The higher states require 
too much excitation energy and belong to different degeneracy multiplets. The initial 
electronic wavefunction is a clear superposition of lscr and 2pcr two-center states. 
One thus expects that only these states are relevant for the calculation. In fact, at 
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these energies the population of the 2p atomic state in charge exchange is much less 
than the population of the Is atomic state. These assumptions are well supported by 



the calculations of Grande and Schwietz |i~0||, who have used a dynamical approach 



based on target-centered wavef unctions. In their approach one has to include a great 
amount of target-centered states in order to represent well the strong distortion of 
the wavefunction as the projectile closes in the target. We also have assumed that 
the proton follows a classical trajectory determined by an impact parameter b. 

Eq. (HD does not look like the usual form of coupled-channels equations in the 
theory of the time-dependent Schrodinger equation. But we can put it in such form 
by rewriting the equation as 

d ( a + \ ( V + + E iW \(a+\ (R . 

where the indices + and — refer to the lscr and 2pcr states, respectively, E = —13.6 
eV, V± (t) = E±(t) - E , and 

_ k {^ + \dV p /dt\^_) = <jw (t)\dV p /dt\y 2pa (t)) 
W{t) ~ n E + (t)-E.(t) ~ n E lsa (t)-E, w (t) ■ {<) 

In this form, the potentials V± (i) and W (i) act like potentials in usual coupled- 



channels equations. We use the formalism of Teller |TT| to calculate the wavefunctions 
\fr± (R) at different inter-proton distances, R(t), corresponding to a particular time 
t. The static Schrodinger equation is solved in elliptical coordinates. This yields 
two coupled differential equations which can be solved by expanding the solutions in 
Taylor series. A set of recurrence relations is obtained for the expansion coefficients 
when the boundary conditions are used. The energies Ei sa (R) and E 2pa (R) are 
obtained by adjusting the constant which separates the two coupled equations to 
its correct matching value. 

When t — > ±00, V± — > and W — > 0. The initial state, an electron localized 
in the target can be written in terms of the degenerate symmetric, \& + = ^i scr , and 
anti-symmetric, = \I / 2po-, states: 

$ T = —= + + *_) , at t — > -00 , (8) 
v2 

where both $t, and are normalized wavefunctions. If the electron is localized 
in the projectile, the wavefunction <3> p = — ^_) /y/2, when t — >■ —00 is used. 
We will consider only the condition of eq. , namely, an electron localized at the 
target at t — > —00. These relations are well known quantum mechanical results; 
the asymptotic two-center wavefunctions can be written as combinations of target- 
and projectile-centered ls-wavefunctions: ^± = ($ p ± <Ev) / 

Starting with a target localized electron we assign the initial conditions a± = 1/ y/2 
at t — > —00 and solve the equation numerically. Although at t — >■ —00 the 
probabilities |a±| 2 remain very close to 1/2, the amplitudes a± acquire phases which 
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change the relative population of the projectile and the target Is state. We correct 
for energy conservation which feeds the increasing binding energy of the electron back 
to an increasing relative motion energy of the two protons as they come closer. This 
is specially important as E p becomes of order of hundreds of eV, and smaller. In 
figure 1 we show the time dependence of V± (t) and W (t) for E p = 10 keV and a 
nearly central collision, b = 0.1 A. One observes that the potentials V± (t) extend 
much farther out than W (t) . Moreover, we find that as E p decreases the potential W 
decreases faster than the projectile's velocity, v p . This is mainly due to the derivative 
of in eq. (|7|). At E p ~ 100 eV the potential W loses its relevance as compared 
to V±, which have no dependence on v p . This becomes clear in figure 2. In this 
figure we show the exchange probability as a function of the impact parameter for 
two projectile energies. The solid line is the full solution of eq. @ . The dashed line 
is the approximation obtained when we set W = in eq. (^) . In the later case, the 
equations decouple and it is straightforward to show that the exchange probability is 
given by 
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^a ± (oo) ($t|*±(oo)) 
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(9) 



At E p = 10 keV there is an appreciable difference between the full calculation and 
the approximation @. But, for E p = 100 eV the results are practically equal, except 
for very small impact parameters at which the potential W still has an effect. 

One observes that the exchange probability is not constant at small impact pa- 
rameters, but oscillates wildly around 0.5, specially for low projectile energies. One 
might naively assume that because the collision is almost adiabatic, the system loses 
memory of to which nucleus the electron is bound after the collision. Thus, for small 
impact parameters one would expect a 50% probability of finding the electron in one 
of the nuclei at t — oo. However, this is not what happens. From eq. @ we see that 
minima of the probability occur for impact parameters satisfying the relation 

/oo 
[£L(t) - E + (t)}dt = 2nh(n + 1/2), n = 0, 1, 2, N. (10) 
-oo 

This relation looks familiar, of course. It simply states that the interference between 
the lscx and the 2p<r states induces oscillations in the exchange probability. The 
electron tunnels back and forth between the projectile and the target during the 
ingoing and the outgoing part of the trajectory. When the interaction time is an exact 
multiple of the oscillation time, a minimum in the exchange probability occurs. The 
average probability over the smaller impact parameters is indeed 0.5. As the impact 
parameter decreases from infinity, the first maximum in the exchange probability 
indicates the beginning of the region of strong exchange probability. One sees that at 
low proton energies this starts at b ~ 3 A. The size of the hydrogen atom is about 0.5 
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A and thus the electron travels in a forbidden region (tunnels) of about 2 A from the 
target to the projectile. This is possible because of the strong interference between 
the lscr and the 2p<j states, which for some trajectories satisfy the quantum relation 

(0)- 

To obtain the stopping power we need the total cross section for charge exchange, 
a = 2tt J P exc hbdb. This is shown in figure 3. The solid line is the full coupled-channels 
calculation, while the dashed line uses approximation (|^) for the exchange probability. 
We observe that the approximation (^) reproduces well the full calculation even at 
the highest energies. The reason is that the potential W is always smaller than V± 
for large impact parameters which have more weight on the integral cross section. 



We also compare our calculations with the lowest energy data of McClure [13|. The 



formalism developed here is inappropriate for energies in the tens of keV range and 
higher, as the projectile velocity becomes comparable to or higher than the electron 
velocity. This implies that two-center states with higher energy and even continuum 
states (ionization) should be included in the calculation. For E p — > 0, the charge 
exchange cross section becomes the constant value a (E p = 0) = 37.88 xlO -16 cm 2 . 
This happens because, when E p — > and as the projectile nears the targets, the 
increasing electron binding in the two-center system acts as a push in the relative 
motion energy to compensate for energy conservation. The average result is that 
the cross section for charge exchange becomes approximately constant for projectile 
energies of tens of eV and below. 

In figure 4 we show the stopping cross section of the proton. The stopping cross 
section is defined as S = J2i A-E* <7j , where AEi is the energy loss of the projectile 
in a process denoted by i. The stopping power, Sp = dE/dx, the energy loss per unit 
length of the target material, is related to the stopping cross section by S = Sp/N, 
where N is the atomic density of the material. In our charge exchange mechanism the 
electron is transferred to the ground state of the projectile and the energy transfer is 
given by AE = m e v p /2, where v p is the projectile velocity. Assuming that there is 
a few free electrons in the material (e.g., in a hydrogen gas) only one more stopping 
mechanism at very low energies should be considered: the nuclear stopping power. 
This is simply the elastic scattering of the projectile off the target nuclei. The pro- 
jectile energy is partially transferred to the recoil energy of the target atom. The 
stopping cross section for this mechanism has been extensively studied by Lindhard 
and collaborators (see, e.g., ref. ]T4|). The nuclear stopping includes the effect of 
the electron screening of the nuclear charges. 

The dotted line in figure 4 gives the energy transfer by means of nuclear stopping, 
while the solid line are our results for the charge-exchange stopping mechanism. The 
data points are from the tabulation of Andersen and Ziegler ||. We see that the 
nuclear stopping dominates at the lowest energies, while the charge-exchange stopping 
is larger for proton energies greater than 200 eV. Since we neglect the difference 
between molecular and atomic hydrogen targets, there is a limitation to compare 
our results with the experimental data. But, the order of magnitude agreement is 
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very good in view of our simplifying assumptions. We do not consider the change of 
the charge state of the protons as they penetrate the target material. The exchange 
mechanism transforms the protons into H atoms. These again interact with the target 
atoms. The can loose their electron again by transfer to the Is state of the target 



The best fit to our calculation for the stopping power for proton energies in the 
range 100 eV - 1 keV yields S ~ t^' 35 . This contrasts with the extrapolation S ~ v p , 
based on the Andersen- Ziegler table. But, this discrepancy is much less than the one 
obtained by Golser and Semrad || for helium targets, who found a stopping power for 
protons S ~ t> 3,34 for protons in the energy range of 4 keV. No data at lower energies 
are available in this case. But, the Golser and Semrad data, for proton energies above 
3 keV, firmly indicate that a high power dependence on the projectile velocity will 
be also valid at lower energies, in contrast the predictions from the Andersen and 
Ziegler tables ||. One cannot extend our calculations to helium targets as the initial 
wavefunction cannot be described in terms of a simple sum of two-center states. A 
much larger two-center basis is necessary. Since the electrons in the helium target 
are more bound than in the proton, the charge-exchange probability must be much 
smaller than in the case of hydrogen targets. One thus should indeed expect a much 
stronger dependence of the stopping on the projectile velocity. At very low energies, 
of the order of some hundreds of eV, the stopping cross section should be entirely 
dominated by nuclear stopping, even more than for hydrogen targets. 

The p + p — > d + e + + v e reaction is a very important one occurring in, e.g., our 
sun. But, it proceeds via the weak interaction and its cross section is extremely small 
for studies under the laboratory conditions 0. Fortunately, a good theoretical 
model exists for this reaction fIB| . Other reactions could be strongly influenced by 
the stopping power of protons and deuterons due to the charge-exchange mechanism. 
They can be relevant for the study of d + D reactions in stellar interiors and fusion 
reactors. Another application is the D(p, 7) 3 i7e reaction which is important for the 
hydrogen burning in stars. In our sun the most effective energy of this reaction is 
E c .m. = 6.5 ± 3.3 keV at T = 15 x 10 6 K. At this energy one expects that the charge- 
exchange stopping cross section should be as important as the ionization cross section. 
Experimental data exist at the lowest energy value of 16 keV [15|, [17]]. Although the 
extrapolation based on theory appears to be under control in this case, it is worthwhile 
to consider a better study of the stopping power for this reaction. The steep rise of 
the fusion cross sections at astrophysical energies amplifies all effects leading to a 
slight modification of the projectile energy [JTSf. Our results show that the stopping 
mechanism does not follow a universal pattern for all systems. This calls for improved 
theoretical studies of charge-exchange effects and for their independent experimental 
verification. 
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Figure Captions 

Fig. 1 - Time dependence of the interaction potentials V± (t) and W (t) for E p — 10 
keV and a nearly central collision, b = 0.1 A. 

Fig. 2 - The exchange probability as a function of the impact parameter for two 
projectile energies. The solid line is the full solution of eq. @ . The dashed line is 
the approximation obtained when we set W = in eq. (0) . 

Fig. 3 - The solid line is the full coupled-channels calculation for the charge- 
exchange cross section, while the dashed line uses approximation (|]) for the exchange 
probability. The experimental data are from McClure [13[ . 



Fig. 4 - The stopping cross section of protons on H-targets. The dotted line in 
gives the energy transfer by means of nuclear stopping, while the solid line are our 
results for the charge-exchange stopping mechanism. The data points are from the 
tabulation of Andersen and Ziegler M. 
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